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We study a system of interacting spinless fermions in one dimension which, in the absence of 
interactions, reduces to the Kitaev chain [A. Yu Kitaev, Phys.-Usp. 44, 131 (2001)]. In the non¬ 
interacting case, a signal of topological order appears as zero-energy modes localized near the edges. 
We show that the exact ground states can be obtained analytically even in the presence of nearest- 
neighbor repulsive interactions when the on-site (chemical) potential is tuned to a particular function 
of the other parameters. As with the non-interacting case, the obtained ground states are two-fold 
degenerate and differ in fermionic parity. We prove the uniqueness of the obtained ground states 
and show that they can be continuously deformed to the ground states of the non-interacting Kitaev 
chain without gap closing. We also demonstrate explicitly that there exists a set of operators each 
of which maps one of the ground states to the other with opposite fermionic parity. These operators 
can be thought of as an interacting generalization of Majorana edge zero modes. 

PACS numbers: 71.10.Pm, 71.10.Fd, 75.10.Pq 


I. INTRODUCTION 

In recent years Majorana zero modes have been in 
the focus of research in condensed matter physics.^i®! 
Experimental signatures of Majorana zero modes have 
been observed in the tunnelling conductance in hv- 
brid superconductor-semiconductor nanowire systems.^^HSl 
However, an unambiguous identification of Majorana 
particles remains open.l^ Very recently an important 
steppe towards such an identification was taken by scan¬ 
ning tunnelling microscopy experiments on ferromagnetic 
atomic chains on a superconductor, which allowed to ob¬ 
tain spatially resolved signatures showing localization of 
the zero modes at the edges as predicted by theory. 

The prototypical toy model possessing Majorana zero 
modes is the Kitaev chain with open boundaries,^ a one¬ 
dimensional tight-binding model for spinless fermions in 
the presence of p-wave superconducting pairing. This 
model possesses a topological phase with two-fold degen¬ 
erate ground states which cannot be distinguished by any 
local order parameter. There exist zero-energy modes 
that commute with the Hamiltonian and anticommute 
with the fermionic parity; these are the Majorana zero 
modes. They are exponentially localized near the bound¬ 
aries. As is well-known (see e.g. Ref. IT^ the Kitaev 
chain can be mapped to the one-dimensional transverse- 
field Ising chain via a non-local Jordan-Wigner transfor¬ 
mation. In the resulting spin variables, the topological 
phase reduces to the ordered phase showing spontaneous 
magnetization while the fermionic parity maps to a string 
of spin operators. 

The recent surge of interest in condensed-matter real¬ 
izations of Majorana fermions includes the investigation 
of the effects of disordeiJi^®^ or dimerizatiorP^ on the 
topological phase, the study of nanostructures possess¬ 


ing Majorana zero modes, strongly correlated 
systems showing Kondo physics Furthermore, gen¬ 
eralizations to higher symmetries including sup ersym me- 
trjIMMI and parafermion modes were analyzed. 

Here we foc us on interaction effects. Motivated 
by the proposaP^HIl to realize systems hosting Ma¬ 
jorana zero modes in semiconducting nanowires and 
the subsequent experimental works,™ 2 l niost previous 
studies of interaction effects have focused on nanowire 
setups. This included the original spinful nanowire 
system with spin-orbit coupling, a Zeeman field and 

proximity to a superconductor,simila r mul ti-band 

nanowires,!^ and effectively spinless systemJSE^ includ¬ 
ing helical wire^^ and two-chain ladders.l^ Further¬ 
more, interaction s directly in the Kitaev chain were 
studied,l^ 2 l 3 lIi 2 lilHIIl .^^Fich were also shown to be imple- 
mentable in an array of Josephson junctions,!^ a realiza¬ 
tion allowing to reach strong interaction strengths as well 
as good control of the parameters. From all these works 
a physical picture emerged revealing two main effects of 
interactions: On the one hand, interactions suppress the 
bulk gap and thus decrease the stability of the topolog¬ 
ical phase, while on the other hand, repulsive interac¬ 
tions broaden the chemical-potential window over which 
Majorana zero modes exist. Which of these two effects 
dominates or is more relevant depends on the specific re¬ 
alization at hand. General aspects of interaction effects 
on topological phases were also analyzed in Refs. [52H581 

In this article, we further investigate the Kitaev chain 
in the presence of repulsive interactions between the spin¬ 
less fermions. Employing the Jordan-Wigner transfor¬ 
mation, this model can be mapped to an XYZ Heisen¬ 
berg chain, all owing to obtain a wealth of results on the 
phase diagram .ESESljjej. 0 ^ however, we are not interested 
in the phase diagram but directly in the two-fold degener- 
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ate ground states in the topological phase. Using classic 
results in the spin chain literaturep^ffiB we find that the 
exact ground states of the model can be obtained when 
the chemical potential is tuned to a particular function of 
the other parameters: the hopping amplitude, the p-wave 
pairing gap, and the strength of the repulsive interaction. 
In this special case, the ground states can be written in 
a simple product form. The ground-state degeneracy is 
a necessary but not sufficient condition for topological 
order and the existence of Majorana zero modes. In or¬ 
der to prove the existence of topological order, we show 
that there exists a smooth path that interpolates between 
the solvable Hamiltonian and the non-interacting Kitaev 
chain in the topological phase, along which the ground 
states remain unchanged. We also prove rigorously that 
the spectral gap above the ground states does not van¬ 
ish along the entire path. Then, thanks to the results of 
Fidkowski and Kitaev,!^ it follows that the topological 
order of the non-interacting Kitaev chain persists along 
the path. To provide further evidence for the topological 
order, we prove the existence of a set of operators each 
of which maps one of the ground states to the other with 
the opposite fermionic parity. The operators obtained 
(i) are hermitian, (ii) anticommute with the fermionic 
parity, and (iii) are localized near the edges. Therefore, 
they can be regarded as an interacting generalization of 
the Majorana zero modes. In fact, in the absence of the 
interactions, they exactly commute with the Hamiltonian 
and reduce to the standard Majorana zero modes in the 
original Kitaev chain. 

The paper is organized as follows. In Sec. [Hj we give a 
precise definition of the model and list the symmetries of 
its Hamiltonian. In Sec. |HI[ we first introduce the notion 
of frustration-free Hamiltonians. Then we show how one 
can find the ground states of a chain of arbitrary length 
from the exact results of the two-site problem. The ex¬ 
pressions for the ground states are also given explicitly. 
In Sec. EYi we discuss the topological order of the present 
system. We present a theorem about the continuous de¬ 
formation between the interacting and non-interacting 
Kitaev chains that share the same ground states. We 
show how the theorem follows from a lemma about the 
auxiliary free-fermion problem. We also present the proof 
of the lemma and derive explicit expressions for the Ma¬ 
jorana zero modes in the auxiliary free-fermion problem. 
Concluding remarks are presented in Sec. [y] In App. 
we summarize the relation between Majorana and com¬ 
plex fermions. The explicit expression for the Hamilto¬ 
nian in terms of Majorana fermions is also given. In 
App.ll we show the mapping between the fermionic 
Hamiltonian and the XYZ chain in a magnetic field. 
We also show how the frustration-free condition and the 
ground states obtained translate into the spin language. 
In App. C, we present a derivation of the inequality used 
in Sec. HI. A few examples of correlation functions in the 
ground states are also provided. In App. D, we extend 
the theorem to the case where the couplings are spatially 
inhomogeneous. In App. we present a detailed expo¬ 


sition of the eigenvalue problem related to the auxiliary 
free-fermion problem. 

II. MODEL 

We consider a system of spinless fermions on a chain 
of length L with open boundaries. For each site j = 
1 , 2, • • • , L, we denote by c| and Cj the creation and the 
annihilation operators respectively. As usual, the number 
operators are defined by rij := c^jCj. 

A. Hamiltonian 

We consider the Hamiltonian of interacting spinless 
fermions described by 

L-l 

H + c]+iCj) + A{cjCj+i + c]+icj) 

i=i 

^ L L-l 

- -1) + U ^(2nj - I)(2n^+1 - 1), 

1=1 i=i 

( 1 ) 

where t is the hopping amplitude and A the p-wave par¬ 
ing gap, which is assumed to be real. Without loss of 
generality, we can assume that t > 0 because the case 
with t < 0 can be achieved by local unitary transfor¬ 
mations: Cj —> —i{—iycj. We can further assume that 
A > 0 because the case with A < 0 can be achieved by 
Cj —> icj for all j. Here, ^j is the on-site (chemical) po¬ 
tential and C/ > 0 is the strength of the nearest-neighbor 
repulsive interaction. It is convenient, for later purposes, 
to keep a site-dependent on-site potential jij . In the ab¬ 
sence of the interaction ([/ = 0), the model reduces to the 
Kitaev chain,E2 in which Majorana edge zero modes oc¬ 
cur provided that it is in the topological phase. Thus one 
can think of the interacting system as Majorana fermions 
with a quartic interaction. An explicit expression for the 
Hamiltonian in terms of Majorana fermions is given in 
App. 0 One can also express the Hamiltonian in terms 
of spin-1/2 operators via a Jordan-Wigner transforma¬ 
tion. The corresponding model turns out to be the XYZ 
spin chain in a magnetic field (see App. for more de¬ 
tails) . 

B. Symmetries 

Let us first consider the symmetries of the Hamiltonian 
H. Because of the presence of the pairing term, H does 
not conserve the total fermion number F := ’^^^I'rij, 
i.e., [H,F] y 0. However, the fermion number mod¬ 
ulo two is conserved, [iL, (—1)'^] = 0. The Hamiltonian 
respects time-reversal symmetry, i.e., is invariant under 
complex conjugation. 
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In addition, when fij = 0 for all j, the Hamiltonian is 
invariant under the charge conjugation operation Cj —>■ 
{—lyCj. More precisely, H commutes with the following 
unitary operator 


P = 


n[c,+(-i)^' 


i=i 



( 2 ) 


We also note that this particular case is integrable, be¬ 
cause it can be mapped to the XYZ spin chain (without 
magnetic fields) which was solved by Baxter!^ For the 
case with non-vanishing /r^-’s, the Hamiltonian H is not 
invariant under the above charge conjugation. However, 
we can at least say that H with and that with 

li8,ve the same spectrum. In the following, for 
simplicity, we assume that > 0 for all j. 


III. EXACT GROUND STATES 

In this section we show that the Hamiltonian H is 
frustration-free when the /r^-’s are tuned to a particular 
function of the other parameters (t. A, U), in which case 
the exact ground states are easy to obtain analytically. 

Consider the case where /ii = /iz, = /r/2 and = /r 
(j = 2,3,..., L — 1), i.e., the on-site potentials at the 
boundary sites are half the bulk ones. In this case, the 
Hamiltonian takes the form 

L-l 

H='£hj, (3) 

1=1 

where 

hj = -t{c]cj+i -f -k A{cjCj+i -k c]+icj-) 

- + rij+i - 1) + U{2nj - I)(2nj+i - I). (4) 


such that Cfejvac) = 0 for any k. Since the local Hamil¬ 
tonian hj commutes with the fermionic parity (—1)^, we 
can deal with even and odd sectors separately. In the 
even sector spanned by | o o) and | • •), hj is expressed 
as 

|oo) !••) 

,(even) _ (U + flj 2 "A \ (5) 

1 -A U -[il2)- 

An unnormalized ground state of and the corre¬ 

sponding eigenvalue are 

= I oo)-Pcot^l ••), (6) 

VA2-t-(M/2)2, (7) 

respectively, where Q = arctan(2A/^) which is assumed 
to be in [0,7r]. Note that is non-degenerate as 

long as either A or /i is non-vanishing. Similarly, in the 
odd sector spanned by | • o) and | o •), hj is expressed as 


|•o)|o•) 



( 8 ) 


An unnormalized ground state of and the corre¬ 

sponding eigenvalue are 

= I *0) + I o*), (9) 

=-((7 + t), (10) 

which is non-degenerate except when t = 0. 

The absolute ground state of hj can be obtained by 
comparing and , which become degenerate if 

I P - A2 

= fi* := 4^tU-k - - -. (II) 


Since [hj,hi^] ^ 0 if \j — k\ = 1, the local Hamiltoni¬ 
ans cannot be diagonalized simultaneously. However, 
there may exist some special values of /i at which the 
projection onto the ground-state subspace Pq satisfies 
hjPo = eqPq for all j, where eg is the smallest eigenvalue 
of hj. In other words, the ground states of H minimize 
each hj independently. In such cases, the Hamiltonian 
H is said to be frustration-free. Many known exactly 
solvable models fall into this category. Examples include 
the AfHeck-Kennedy-Lieb-Tasaki model,!^ the Rokhsar- 
Kivelson model,!^ and the Kitaev toric code model.^ 
To search for the condition under which the Hamil¬ 
tonian ([^ is frustration-free, we shall first consider the 
two-site problem, i.e., deriving the spectrum and the 
eigenstates of hj. Let us ignore, for the moment, the 
region outside the bond {j,j -|- 1). Then the Hilbert 
space for the problem is spanned by the four states: 
I o o) := |vac), 1*0) := c]|vac), | o •) := c]+i|vac), 
and I •• ) := c|4_|_]^|vac), where |vac) is the vacuum state 


Note that though = —fj,* is also allowed, we have cho¬ 
sen the p ositive /i for simplicity (see the discussion in 
Sec. HB). At /r = /i*, any linear combination of 
and ) is also a ground state of hj. If we take the 

following particular combinations, we see that they are 
disentangled 


=(1 ± q ; c ])(1 ± ac]_,_;^)|vac) 

= exp(±ac]) exp(±ac]^J|vac), (12) 

where 0* = arctan(2A//r*) and a = •\/cot(0*/2). We 
recall that 9 € [0, tt] has been assumed. In the entire 
Hilbert space, the ground states of hj at fi = fj.* are 
highly degenerate, because any state of the form 

f{c\, ..., c]_i) 5(c]-e2: cl)\vac) (13) 
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is a ground state of hj , where / and g are arbitrary poly¬ 
nomials in c \,and cj+ 2 ’ '-l; respectively. 

In this way, we find that the ground states of hj at /r = 
/i* can be factorized into the product of two states. This 
observation hints at the possibility that we may obtain 
exact ground states of H with g = fi* for any L. In fact, 
H is frustration-free at g = g*. To see this, consider the 
states of the form 

... e±«4|vac). (14) 


with jk — ji = i — 1. Noting that {Oo,(—1)'^} = 0 
and using the relations (—1)^ 

(—1)-^ finds 

(Oo)even = (Oo)odd = 0, (19) 

for any Oq. On the other hand, for Oq, one finds that 
the difference between (Oe)even and (Oe)odd is bounded 
from above as follows: 

I (Oe)even ~ (Oe)odd| ^ C||Oe||e ( 20 ) 


Here the coefficient has been introduced so that 
(d'p^^= 1. Since exp(±acj(.) commutes with hj 
unless k ^ j^j + I, one can easily verify that min¬ 

imizes each hj independently at g = g*. Therefore, H is 
frustration-free at a = g* and are its exact ground 

states with energy^ 

Eo = -iL-l){U + t). (15) 


where ||Oe|| denotes the operator norm of Oe, 




1 ,, i-pcot(r/2) 

—;—r with n = -^, 

Injiyl 1 — cot(0*/2) 


and the coefficient C is given by 


O l + l^l 

l-l/r;2 


( 21 ) 


( 22 ) 


A similar factorization of the ground states has been 
found in spin-1/2 chains in a magnetic fieldfact, 
the condition under which the Hamiltonian is frustration- 
free is the same as the exactly solvable condition of the 
XYZ spin chain in a magnetic field (see App.j^for more 
details). 

It should be noted that are not orthogonal 

and are not eigenstates of the fermionic parity (—1)'^, 
because each contains both even and odd num¬ 

bers of fermions. The ground states that are eigen¬ 
states of (—1)'^ can be obtained by taking appropriate 
linear combinations of |4 'q^^). Let us decompose the en¬ 
tire Hilbert space H as T-L = "Hi®™") 0 where 

(even/odd) QQnsists of the states with even/odd numbers 
of fermions. Then, noting that [F, c]] = c] and hence 
(—1)^ exp(±ac]) = exp(0acj)(—1)^, we find that the 
following states 


Clearly, the result shows that the states and 

cannot be distinguished by any local measure¬ 
ment for large L. A detailed proof of the inequality Eq. 
^(|20[ ) as well as some explicit examples are given in App. 

Let us finally comment on the location of the solvable 
case Eq. (Ill in the phase diagram. The phase diagram 
of the model described by the Hamiltonian Eq. ([^ with 
A = t an d gj = g (4 = 1,2, ...,L) has been worked out 
previousljB^® Fig. [l] shows the phase diagram for this 
particular case. Since the spectrum of the Hamiltonian 
is invariant under sending g — >■ —/i, we only show the 
region where g/t > 0. The origin {U/t,g/t) = (0,0) 
corresponds to the original non-interacting Kitaev chain, 
which is in a topological phase. When A = t, the 
frustration-free condition reads 


g = 4^/U^+tU, 


(23) 


|^(even)) ^ ^ ^ (16) 

|vI;(°dd))^|^(+)^_|^(-))g^(odd), ( 17 ) 

are orthogonal to each other. As we will prove in the next 
section, there is no other ground state of H with g = g* 
except for these two states. Furthermore, for large L, 
they cannot be distinguished by any local measurement. 
To illustrate this, let us consider the expectation values 
of local operators in For notational conve¬ 

nience, we write 


(^(par)|^(par)) 


(18) 


where par is either even or odd. Let Oq/Oq be a lo¬ 
cal operator consisting of even/odd number of creation 
and annihilation operators. We assume that Oq/Oq is 
supported on a set of lattice sites ji < j 2 < ■ ■ ■ < jk 


which is shown by the red (dashed) line in Fig. Clearly, 
this line is in the topological phase. We note that both 
the frustration-free line and the phase boundary between 
the trivial and topological phases approach g = iU in the 
infinite-{7 limit. In this limit, the model can be mapped 
to a model of hard-core spinless fermions, which is solv¬ 
able using the free-fermion method.l^ 


IV. TOPOLOGICAL ORDER 

In this section we discuss the topological order of the 
ground states obtained in the previous section. We ex¬ 
plicitly construct a one-parameter family of Hamiltonians 
that interpolates between interacting and non-interacting 
Kitaev chains, the latter of which exhibits topological or¬ 
der. The spectral gap above the ground state and the ex¬ 
istence of Majorana zero modes that are localized around 
the edges are proven rigorously. 
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FIG. 1: (color online). Phase diagram of the interacting Ki- 
taev chain Eq. Q with A — t and /ij = fi. The region 
fi/t < 0 is not shown since the phase diagram is invariant 
under sending jj, ^ —jj,. The red (dashed) line indicates the 
frustration-free line Eq. (231, along which the exact ground 
states can be obtained in closed form. ICDW stands for the 
phase where an incommensurate charge-density wave state is 
a ground state. 


A. Smooth path connecting interacting and 
non-interacting Hamiltonians 


Let us introduce the Hamiltonian that interpolates 
continuously between the interacting and non-interacting 
model. As discussed in the previous section, the ground 
states depend only on 0* = arctan(2A//r*) at 

^ = 11 * (frustration-free case). Therefore, they are also 
ground states of the following one-parameter family of 
the Hamiltonians 

L-l 

H{s) = Y. hj{s) with s > 0, (24) 

i=i 

where the local Hamiltonians are given by 


hj{s) = -{CjCj+i -b h.c.) -b (1 -b s) sm9*{cjCj+i + h.c.) 
— (1 -b s) cos 9* {rij + Uj+i — 1) 

+ ~ ~ 2) ■ 

Here we have set t = 1. The last constant was added 
so that the ground-state energy is equal to zero. Since 
we have assumed non-negative A and /r*, 0 is in [0,7r/2]. 
When /r = /i*, the original Hamiltonian H in Eq. © is 
the same as H{2U) up to a trivial constant, which can be 
verified by direct calculation. On the other hand, since 
the second last term in Eq. (25) vanishes when s = 0, 


iL(0) reduces to the non-interacting Kitaev Hamiltonian 
with pairing A/(l-b2t7) and on-site potential /i*/(l-b2[/) 
which we denote by Hq in the following. Thus we see that 
the Hamiltonians H and Hq are smoothly connected to 
each other. Along the entire path 0 < s < 2U, the ground 
states remain unchanged.!^ This suggests that the two 
Hamiltonians H and Hq are in the same topological phase 
as long as the energy gap does not close along the path. 

Now we show the existence of an energy gap of H{s). 
To this end, we prove the following 


Theorem: For all s > 0 and 0* 7^ 0, (i) the ground state 
of the Hamiltonian H{s) is unique up to double degener¬ 
acy, (ii) H(s) has a uniform (independent of the length 
of the chain) spectral gap above the ground state. 

The proof of the theorem relies on the following 

Lemma: For 9* 7^ 0, (i) the ground state of the non¬ 
interacting Hamiltonian H(0) is unique up to double de¬ 
generacy, (ii) H{0) has a uniform (independent of the 
length of the chain) spectral gap above the ground state. 

This lemma can be proven by directly computing the 
spectrum of H{0), which will be given in the next sub¬ 
section. We shall now show how the theorem follows from 
the lemma. 

Proof of Theorem: It is convenient to introduce the fol¬ 
lowing operator 


9 


Qj = T cos ^ {cj + Cj+i)(c] - cj-+i)(c] -b cj-+i) 


A n A 


9* 


+ ^ sin y (cj - cJ_^i)(cj -b Cj+i){cj - Cj+i), (26) 

which is a sum of triple products of fermion operators.!^ 
In terms of Qj, hj{s) in Eq. (25) is rewritten as 


“ QjQ^j T (1 T s)Q\Qj 




(27) 


which can be verified by a tedious but straightforward 
calculation. The above form of hj{s) is manifestly pos¬ 
itive semidefinite for all s > — 1. It is then clear that 
hj{s) > hj(0) for s > 0. Here, we write A > B 
to denote that A — B is positive semidefinite.Since 
the inequalities for the local Hamiltonians hold for all 
j = 1,2,..., L — 1, we arrive at the inequality 


H{s) > H{0), for s > 0. (28) 


The states are annihilated by both Qj and Q] for 

all j, and hence are ground states of H{s). 

We shall next compare the energy eigenvalues of H(s) 
and H{0). Let En{s) be the nth eigenvalue of H{s). Here 
the energy eigenvalues are arranged in non-decreasing or¬ 
der: 0 = A’i(s) = ^2(5) < E^{s) < ... < En{s) < .... 
Note that H{s) has at least two zero-energy ground 
states. It follows from min-max principl^^ that the in¬ 
equality Eq. (28) implies E„(s) > A„(0) for all n. There¬ 
fore, the spectral gap is a non-decreasing function of s, 
and the eigenstates of H{s) continuously connected to 
those of H{0) with positive energy never join the zero- 
energy manifold with increasing s. This, together with 
the lemma, proves the theorem. ■ 


Several comments are in order. Our theorem implies 
that the interacting Hamiltonian H{2U) and the non¬ 
interacting one H{0) are adiabatically connected to each 
other without gap closing. Thanks to the work of Fid- 
kowski and Kitae\El, this suffices to show that H{2U) 
is topologically non-trivial, because the non-interacting 
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H{0) is in a topological phase. However, some care is 
needed here. As shown in Refs. [5353], the classifi¬ 
cation of free-fermion topological phases breaks down 
in one dimensions in the presence of interactions. The 
non-interacting Kitaev chain with time-reversal symme¬ 
try belongs to the BDI symmetry class, different phases of 
which are characterized by Z integers. This Z classifica¬ 
tion is broken down to Zg in the presence of interactions. 
One might think that there is a possibility that a topo¬ 
logically non-trivial phase is adiabatically connected to 
the trivial phase with no Majorana edge zero modes via 
interactions. This is, however, not the case here, because 
the topological index is basically the number of Majo¬ 
rana edge modes localized at one end of the chain. For 
H{0), the number of Majorana edge zero modes at each 
edge is 1, which is apparently different from 0 modulo 8. 
Therefore, the one-parameter family of H (s) with s > 0 is 
in a topologically non-trivial phase which is adiabatically 
connected to the one in the non-interacting classification. 

We also comment on the stability of the spectral gap 
against small perturbations. Let us first consider the 
effect of the inhomogeneity arising from the on-site po¬ 
tential term. As noted in Sec. m the on-site poten¬ 
tials at the boundary sites are half the bulk one in the 
Hamiltonian Eq. (|^. To see whether the spectral gap is 
robust against the perturbation that restores the homo¬ 
geneity to some extent, consider the perturbed Hamilto¬ 
nian H' = H + V with ^ = fi* and 

E =-yCni-f rii), (29) 

where 6 fi is assumed to be positive. Let En{2U,Sfi) be 
the nth eigenvalue of H'. The condition for the existence 
of the spectral gap reads E^{2U, Sfi) > £’2(217, 6 fj,). From 
Weyl’s theorem (see Theorem 4.3.1 in Ref. EH), one finds 
£„(2£) - 6(1 < E„{2U,S(i) < £„(2t7) for all n. This 
implies that iJ' has a spectral gap if Sfi < £3(217). Since 
our Theorem ensures that £3(217) is strictly positive, we 
conclude that the spectral gap is robust against V for 
small enough 6 / 1 . A more quantitative estimate is possi¬ 
ble, because the lower bound for £3(217) is obtained as 
E 3 I 2 U) > £3(0) = 2(1 — COS0*) in the next subsection. 
Using this value, we find that the spectral gap of H' is 
non-zero provided 


6(1 <2 


Vi 


1-^217)2-A2 
1 + 2U 


(30) 


A more interesting question is whether the spectral gap 
of H{s) is robust against small but global perturbations. 
For quantum spin systems, rigorous perturbation theo¬ 
ries for th e sta bility of gapped ground states have been 
developed! ^2 | 73 | xhough the existing methods usually as¬ 
sume translation invariant Hamiltonians and do not im¬ 
mediately apply to the present system, we expect that 
some modification of them is likely to apply particularly 
to the region around the point (A, 17) = (0,0) corre¬ 
sponding to the classical Ising chain. A thorough anal¬ 


ysis is, however, beyond the scope of the present study 
and is left for future work. 

Finally, we comment on the generalization of our model 
to the case where the couplings are allowed to vary spa¬ 
tially. One can, in fact, construct such a model without 
changing the ground states of the Hamiltonian H(s). The 
key to the construction of such a model is to note that 
a sum of operators each of which is a polynomial in Qj 

and Q] (j = 1, 2, — 1) annihilates the states 

If this sum is positive semi-definite, then these states are 
the lowest eigenstates of the sum with zero eigenvalue. 
One such example is the following Hamiltonian: 

L-l 

= E (31) 

i=i 

where aj and /3j are arbitrary positive numbers. Each 
local interaction can be either repulsive or attractive, de¬ 
pending on I3j/aj. More precisely, it is repulsive (at¬ 
tractive) if Pjlaj > 1 (0 < Pj/aj < 1). For the above 
Hamiltonian, we can prove a result similar to the theo¬ 
rem above, and hence the presence of topological order. 
A precise statement and its proof are provided in Ap¬ 
pendix]^ We emphasize that the above example demon¬ 
strates that the couplings in the Hamiltonian need not 
necessarily be uniform to obtain topological order. 


B. Proof of Lemma 


Let us prove the lemma by explicitly computing the 
spectrum of H (0). Though i7(0) is quadratic in c and , 
it is not quite straightforward to diagonalize it because of 
the boundaries. Remarkably, however, we find that the 
entire spectrum of H (0) can be obtained analytically. To 
show this, we first rewrite i7(0) in terms of Majorana 
operators Oj = Cj +Cj and bj = {cj — c\)/\, the properties 
of which are summarized in App. The Hamiltonian 
reads 


£( 0 ) 


L 

7) ^ ) Bj^kajbk -|- (£ — 1 ), 

j,k=i 


(32) 


where the real matrix B is given by 


/ c 1-s \ 

1 H" s 2c 1 — 5 


B = - 


V 


1 “h S 2c 1 — 5 
1 -bs c / 


(33) 


with s = sin0*, c = cos^"*. Here, matrix elements 
which are zero are left empty. Since the matrix B is 
non-hermitian, it may not be diagonalizable, but can be 
written in the singular-value decomposition (SVD) form 
B = UAV'^ with A = diag(ei, £2,..., ez,), where orthog¬ 
onal matrices U and V diagonalize BB^ and B^B, re¬ 
spectively. The singular values e„ (n = 1, 2,..., £) can be 
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obtained as non-negative square roots of the eigenvalues 
of BB'^. 

To see the relation between these singular values and 
single-particle energies of H{0), we introduce the follow¬ 
ing new set of Majorana operators 

L L 

~ (34) 

i=i i=i 

which satisfy (a')] = a', (5')] = 6), {o', a'^.} = {6', foj,} = 
2Sj^k, and {a',6(,} = 0 for all j,k = 1,2, It then 

follows from the SVD of B that H{Q) is rewritten as 

L 

i7(0) = ^eX&; + (i-l) (35) 

71 — 1 

= (36) 

r> —1 ^ 


where /„ := (a(j + ih'^)/2 are complex fermions satisfying 
{fm, fl} = 5m,n- Therefore, the set of singular values of 
B is exactly the same as the single-particle spectrum of 
H{Q). To get e„, we need to diagonalize BB^. Since this 
matrix is pentadiagonal, a direct diagonalization by ana¬ 
lytical means is not feasible. In the present case, however, 
something special happens. In fact, we have 

BB'^ = C'^ (37) 

where C is a real symmetric and tridiagonal matrix given 
by 


/I 


— s 
c 


c 

2 


C = 


V 


\ 

c 

c 2 c 
C 1 + 5 / 


(38) 


A similar factorization was found by Truong and Peschel 
in the study of the Hamiltonian limit of the corner trans¬ 
fer matrix of the same model.^ Because of the presence 
of the boundary terms, the diagonalization of C is still 
non-trivial, but can be done analytically. The details are 
given in App.|^ The exact eigenvalues of C are as follows 



is positive semi-definite, a state annihilated by /„ for all 
n = 1, 2,..., L — 1 is a many-body ground state w ith zero 
energy. The states I'I'q®™”^) and in Eqs. (16 [^) 

are such states, because we have already shown that they 
are zero-energy states of 77(0). Then by acting with /I 
(n = 1,2,..., L— 1) on these two states, one can construct 
many-body eigenstates of 77(0) as follows: 


(40) 

(/i^rK/ln (41) 


where rrij = 0 or 1 for all j = 1, 2,..., L — 1. Note that 
the energy of the above two states is 77 = 
which implies that every energy level is at least two-fold 
degenerate. In this way, we have obtained 2^ number of 
energy eigenstates forming a complete set of orthogonal 
states, which can be verified by noting = 

0. Because of the positivity of e„ for n = 1, 2,..., L — 1, we 
find that and |'I'o°'^'^) are the only zero-energy 

ground states of 77(0). This proves the uniqueness part 
of the lemma. We now turn to prove the existence of 
the spectral gap. The lowest-energy excited states are 
obtained by filling the lowest positive e„. More explicitly, 
they are given by /i_iand We 

thus obtain the spectral gap as 


773(0) = iL-l = 2 


1 — c cos 



(42) 


which is a decreasing function in L when L > 0. In the 
infinite-L limit, we have 773(0) —)■ 2(1 — cos0*) which is 
strictly positive unless 9* = 0. This completes the proof 
of the lemma. ■ 

There is an alternative proof of the uniqueness of the 
two-fold degenerate gr oun d states using the spin- chain 
Hamiltonian (see Eq. ( |B4[ ) in App. H- Erom Eq. ( |B5| ), 
one finds that all the off-diagonal elements of the Hamil¬ 
tonian in the standard Ising basis are non-positive when 
A > 0. Then, noting that the space of eigenstates is sep¬ 
arated into two disconnected sectors with opposite par¬ 
ities cTj), we find that in each sector the Hamil¬ 

tonian satisfies the connectivity condition when A > 0, 
i.e., 9* ^ 0,7r. It then follows from the Perron-Erobenius 
theoreirPS that the ground state of each sector is unique. 
Therefore, the absolute ground state is at most two-fold 
degenerate. But since we have already found two ground 
states, this proves that the ground-state degeneracy is 
exactly two. 


The presence of the zero-energy mode clearly signals the 
topological order in 77(0). The system indeed possesses 
Majorana edge zero modes as we will see in the next 
subsection. 

Let us now prove the uniqueness of the ground state 
and the existence of the spectral gap, which does not 
vanish as L —>■ 00. From the fact that trC = 2{L — 
1), one finds that = 2{L — 1), which together 

with CL = 0 yields 77(0) = Cn/I/n- Since 77(0) 


C. Majorana edge zero modes 

The presence of a zero-energy single-particle state is 
usually a signal of topological order, but does not neces¬ 
sarily mean the existence of Majorana zero modes that 
are localized at the boundaries. Here we show that, for 
77(0), the explicit expressions for the Majorana operators 
can be obtained analytically. Let us first define Majorana 









edge zero modes^^ A Majorana edge zero mode T is an 
operator such that 

• rt =r 

. [H{0),T]=0 

. {(-i)^,r} = o 

• localized near the edge and normalizable as = 1 
even in the infinite-L limit. 


We next show that one can construct such edge modes 
from the left and right null vectors of B in Eq. (33). To 


see this let u = {ui,U 2 , ■■■, ul) be a left null vector of B. 
An elementary calculation shows that 


u = - ■ 


„L-1 


) with r = — 


1+s’ 


(43) 


satisfies uB = (0, 0, • • • ,0). Then from the anticommu¬ 
tation relations of a and b operators, one finds that the 
following operator commutes with H(0) 


L 

rL=Af^r^-^a„ (44) 

where the normalization factor Af has been introduced 
so that (Fl)^ = 1. Similarly, from the right null vector 
of B, one finds that 

L 

(45) 

i=i 

commutes with H{0). Since both Fl and Fr are lin¬ 
ear in a and b operators, they are manifestly hermitian 
and anticommute with the fermionic parity (—1)^. They 
are localized near the edges and normalizable, because 
their amplitudes decay exponentially with distance from 
the boundary unless 9* = 0. Therefore, the operators 
Fl and Fr satisfy all the criteria, and hence are Majo¬ 
rana edge zero modes. A linear combination of Fl and 
Fr gives a non-local complex fermion which is nothing 
but /l = (Fl -l-iFR)/2, corresponding to the zero-energy 
state of the single-particle Hamiltonian. Then from the 
properties of Majorana edge zero modes, together with 
the uniqueness of the ground states of H{0), it follows 
that either Fl or Fr maps one ground state to the other 
with different fermionic parity. 

A special feature of H{0) that is non-interacting and 
frustration-free is that the obtained Majorana edge zero 
modes exactly commute with H(0) even in finite-size sys¬ 
tems. This nice property no longer holds for H{s) with 
s > 0. Nevertheless, they still map one of the ground 
states of H{s) to the other since the ground states re¬ 
main unchanged with the introduction of s. Thus, they 
can be regarded as an interacting generalization of Majo¬ 
rana edge zero modes. This is one of the main results of 
our paper. We note that one can in principle construct 
an operator that exactly commutes with H{s) from the 


adiabatic continuation of Fl or FrI ^^ I ^^ I Though this re¬ 
quires a full diagonalization of H{s) and is not feasible, 
we expect that the resulting operator which is no longer 
linear in a and b has a significant overlap with Fl or Fr. 
We leave a quantitative analysis for future work. 


V. CONCLUDING REMARKS 

In this paper, we have studied the one-dimensional Ki- 
taev chain with nearest-neighbor repulsive interactions. 
We have shown that the Hamiltonian of the model is 
frustration-free when the on-site (chemical) potential is 
tuned to a particular function of the other parameters. 
Under this condition, the exact ground states can be 
obtained in closed form. We have also introduced a 
smooth path between the interacting and non-interacting 
Kitaev chains, along which the ground states remain un¬ 
changed. We proved the following theorem about this 
one-parameter family of Hamiltonians along this path: 
(i) the ground state is unique up to double degeneracy, 
and (ii) there exists a uniform (independent of the chain 
length) spectral gap above the ground state. The the¬ 
orem implies that the interacting Kitaev chain in the 
frustration-free case is topologically non-trivial, because 
it is adiabatically connected to the non-interacting Ki¬ 
taev chain in the topological phase. The stability of the 
spectral gap against perturbations at the boundary sites 
was also discussed. Furthermore, we have demonstrated 
explicitly that there exists a set of Majorana edge zero 
modes each of which maps one of the ground states to 
the other with the opposite fermionic parity. It would be 
interesting to explore the implications of this interacting 
generalization of Majorana edge zero modes in transport 
properties such as the tunneling conductance. Though 
the dynamical Green functions cannot be obtained ana¬ 
lytically even for the frustration-free case, one might be 
able to develop an efficient computational method us¬ 
ing the information of the exact ground states. It would 
also be interesting to see how the results obtained trans¬ 
late into the language of continuum field theories. As 
is well known, an appropriate scaling limit of the XYZ 
spin chain is described by the sine-Gordon quantum field 
theory.l^ We thus expect that the continuum limit of 
the interacting Kitaev chain is also described by a sine- 
Gordon type model with boundaries. We also expect that 
boundary bound state^^ arising in integrable quantum 
field theories are natural candidates for the interacting 
generalization of Majorana edge zero modes in the con¬ 
tinuum limit. The precise correspondence is, however, 
beyond the scope of the present study and is left for fu¬ 
ture work. 
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Appendix A: Majorana Hamiltonian 

For each site j, we define the Majorana fermions by 
aj := Cj -b c], bj := {cj - c])/i. (Al) 

One can easily see that they satisfy the defining relations: 

for all j, A: = 1, 2, L. When written in terms of aj and 
bj, the Hamiltonian H in Eq. Q becomes 

. L-l 

i=i 

. L L-l 

— - b^jO'jbj — U ajbjttj+ibj+i + const. (A2) 
^ j=i 

Clearly, the last term represents a quartic interaction 


Appendix B: Spin-chain Hamiltonian and the 
ground states 


In this appendix we discuss the relation between 
the Hamiltonian H in Eq. and the spin- 

chain Hamiltonian. From the standard Jordan-Wigner 
transformation,!^ 


Cl 


erf+ia 


y 

1 


2 



aj + ia| 


2 


for 2 < j < L, 


we obtain explicit expressions of Majorana fermions in 
terms of spin operators as 


Oi = (Ti , 

bi = al, 



for 2 < j < L, 



for 2 < j < L, 


(HI) 

(B2) 


where (crJ,CTj,cr|) are the Pauli matrices at site j. We 
note here that the fermionic parity in spin variables is 
written as 


{-lf = l[al. (B3) 


Substituting Eqs. (Bl B2) into Eq. (A2|, we have 


L-l 

H = 

i=i 

L 

— ^ Bja^ + const. 
i=i 


- Jyalal^-^ + J^cr^crl+i) 




(B4) 


where the parameters are given by 

L = L = ^. L = C, (B5) 

Recent theoretical proposals to realize spin-1/2 XYZ 
chain in a magnetic held using optical lattice systems 
can be found in Refs. I5DH5^ 

In the spin-chain language, the frustration-free condi¬ 
tion Eq. (Ill reads Bi = Bl = B*/2 and Bj = B* 
{j = 2,3, ..fi- I), where 


B* = -2^{J, + J^){J, + Jy). 


(B6) 


This is exactly what has been found in Refs. [Ml and fHTl 
When written in terms of spin states, the ground states 
Eq. (141 become simple direct product states: 




(B7) 


1=1 


which can be distinguished by the expectation value of 
the local operator crj. We note, however, that crj for 
general j is non-local in the fermions, because it involves 
a string of fermion operators. The ground states with 
hxed (fermionic) parities, i.e., Eqs. (16 17), are no longer 
direct products, but can still be written in the form of 
matrix product states. Their explicit expressions can be 
found in Ref. |M1 


Appendix C: Expectation values of Oe 


In this appendix we present a derivation of the inequal¬ 
ity Eq. ( [^ . We also provide explicit expressions for sev¬ 
eral kinds of correlation functions in the inhnite-L limit. 


1 . 


Derivation of Eq. 


( 20 ) 


Let us hrst prove the inequality Eq. (20). Noting that 
[Oej = 0 and using the relation (—l)^|'kg'^^) = 

we find 


(Cl) 

= (C2) 
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This yields 


2. Correlation functions 


(^e)even — 
(^e)odd — 


l + l/r]^ 

(vl/W|Oe|^ ^+^)-(^[,+Ve|^^^) 

1 — 1/77'^ 


(C3) 

(C4) 


where we have used the fact that |4'o^^) = 1 and 
= \/r]^ (see Eq. (21) for the definition of 
rf). Then we find 


|(f^e)even (^e)odd| 

(l/l^l^)l|Oe||+ (d/^IOel'I^^^) 


<2 


1 — l/rf^ 

where the following inequality has been used 

< llOell- 


(C5) 


(C6) 


As we have shown above, any local Oq has the same ex¬ 
pectation value in and in the infinite-L 

limit; for completeness we finally give explicit expressions 
for some of correlation functions. For density operators 
rij, we have 


^l^(n,)eve„/odd - ^ ^ tan(0V2) ’ 
lim (j^jT^k^eYen/odd ~ 

L^oo 


1 -I- tan(0*/2) 


(C12) 

(C13) 


Similarly one can compute equal-time Green functions. 
For j < k with \j — A:| < oo, we have 


lim (CjC^)even/odd 

L^oo 

sin 9* 


2(1 -I- sin0*) 


1 - cot(6»V2) 
1 -f cot(6»*/2) 


1 k-j-l 


(C14) 


In the following, we focus on 0 < 0* < tt, which implies 


/el 4'. 


(-)^ 
0 / 


, we recall that Oq 


\t]\ > 1. To evaluate (4'g''’^|Oe 
is a local operator supported on a hnite number of sites. 
Since cj and c] commute with Oe unless j = 
one finds 




(C7) 


where 




(l + a2)C2 


Jk 




are the ground states restricted to the i consecutive lat¬ 
tice sites j = ji, ji -|- 1, Then using the Schwartz 

inequality, one finds 






< IlOe 


(C9) 


Substituting this into Eq. (C5), one finds 


which decays exponentially in the distance \ j — k\. 


Appendix D: Inhomogeneous model 

In this appendix we show t hat the model described 
by the Hamiltonian (Eq. ( [^ ) exhibits topological 
order. To this end, we prove the following 

Proposition: For all aj > 0, j3j > 0 (j = — 

Ij and 9* ^ 0, (i) the ground state of the Hamiltonian 
is unique up to double degeneracy, (ii) has a 

uniform (independent of the length of the chain) spectral 
gap above the ground state. 

This can be proven along the same lines as Theorem in 
the main text. 

Proof of Proposition: We first introduce the minimum 
value of the inhomogeneous parameters as 

7 = min U (Dl) 

which is, by definition, strictly positive (7 > 0). Then 
from the fact that QjQ^j > 0 and Q^jQj > 0, we find 


|(Oe)eve„ - (Oe)odd| < 2 ||Og| 


1 — l/rf^ 




Iryl^ 1 — 1/772 
where the second line follows from 

— ^ — f- < —^— for 0 < a; < 1. 
1 — a:^ 1 — a; 


(CIO) 


(Cll) 


Then it is obvious that the desired inequality Eq. (20) 
follows from Eq. (CIO). 


L-l 

> iY.{Q,Q] + Q]Qj) = lH{0), (D2) 

1=1 

where H{0) is the Hamiltonian of the non-interacting Ki- 
taev chain. The Hamiltonians H{0) and share the 
same ground states because they are annihilated 

by both Qj and Q] for all j. 

Let us next compare the energy eigenvalues of and 
H{0). Let be the Tith eigenvalue of Here the 

energy eigenvalues are arranged in non-decreasing order: 
0 = < ■■■ < Elf'^ < .... It then 
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follo ws from min-max principl^^ that the inequality Eq. 
(D2) implies > 7ii'„(0) for all n. This immediately 
implies that E™'' > 0 since we have already shown that 
2^3(0) is strictly positive unless 0*^0 (see Eq. ([4^). 
This completes the proof of property (ii). Property (i) 
also follows from min-max principle and Lemma in Sec. 
IV A. ■ 


Let us next show the existence of topological order 
in the system described by We can construct the 

Hamiltonian that interpolates between and H (0) as 
follows: 


= a;iL“‘'-b(l-a:)7iL(0), (D3) 


{vi{q),V 2 {q), ■■■,VL{q)) be an eigenvector of C. For each 
component, we choose the following ansatz 

u^(<7) =ae‘'?^+/3e-*«, (El) 

where the coefficients a, /3, and the “wavenumber” q 
will be determined from the matching conditions at the 
boundaries. The eigenvalue equation Cv{q) = e{q)v{q) 
can be written in components as 

(1 - s)vi{q) + cv 2 {q) = e{q)vi{q), (E2) 

cvj-i{q) + 2vj{q) + cvj+i{q) = Vj{q), I < j < L, (E3) 
CVL-l{,q) + {1+5)vL{.q) = ^{q)vL- (E4) 


where 0 < x < 1. Using the inequality Eq. (D2), one can 
check that > jff(O) for 0 < a: < 1. Therefore, 

by repeating the proof of Proposition, one finds that the 
spectral gap of Il‘^^(x) does not close along the path 
connecting jH{Q) and This, together with the re¬ 

sult of Fidkowski and Kitae\E21, implies the existence of 
topological order in the system. 

Let us remark that a local interaction between the 
neighboring sites can be either repulsive or attractive in 
depending on Pj/oij. To see this, we rewrite the 
Hamiltonian as 


= E '■i I-1 


i=i 


(D4) 


where hj{s) is given by Eq. (25). It is then clear that 


the local interaction between rij and Uj+i is repulsive if 
Pj/aj > 1, while it is attractive if 0 < (ijjaj < 1. The 
interactions can be made purely attractive by taking a 
set of aj and Pj that satisfy 0 < Pjiaj < 1 for all j. For 
instance, if the parameters are chosen as 


Substituting Eq. (El) into Eq. (E3), one finds that the 


eigenvalue must be of the form: 

e{q) = 2 + 2c cos q. 


(E5) 


Then, from the matching conditions at the boundaries, 
i.e., Eqs. (E2 E4), the consistency condition for q is 
obtained as 


-c-(l-fs)e‘« -c-(l-bs)e-‘« AM _ /O 

(s - 1 - ce‘9)e‘9'^ (s - 1 - ce-'9)e-‘'?^ j [p)~[o 


The above equation has a non-trivial solution if the de¬ 
terminant of the matrix in the left-hand side vanishes. 
This condition reads 


— 4ic(l -|- c cos q) sm{qL) = 0. 


The real wavenumbers are then obtained as 


9=—, n = 1,2,....,L-1. 


(E6) 


(E7) 


a, = 1, pj = l + 2U (j = l,2,...,L-l), (D5) 

the Hamiltonian reduces to H{2U), which describes 
the Kitaev chain with attractive interactions when U < 
0. Note that the condition U > —1/2 is required so that 
H{2U) is positive semidefinite. 


Note that q 
result in v{q) 
solution 


0 and TT are not all owed because they 
0. The equation (E6) also has a complex 


9 = 9o = TT -I- i In 



(E8) 


Appendix E: Eigenvalues of the matrix C 

In this appendix we present a detailed exposition of 
the calculation of the eigenvalues of C in Eq. ( [M| ). We 
first make an ansatz for the eigenvectors. Let v{q) = 


which is a solution of cos go = On® might think 

that the set of solutions obtained is overcomplete because 
q = —go is a solution as well. However, this is not the case 
since both go and —go yield the same v{q). We finally 
obtain the desired eigen valu es Eq. ( [^ by substituting 
the obtained g into Eq. (E5). 
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